Use and Abuse of a Fractional Fokker-Planck Dynamics for 

Time-Dependent Driving 

E. Heinsalu, 1 ' 2 M. Patriarca, 1 ' 2 I. Goychuk, 2 and P. Hanggi 2 

1 Institute of Theoretical Physics, Tartu University, Tahe 4, 51010 Tartu, Estonia 
2 Institut fur Physik, Universitat Augsburg, 
Universitdtsstr. 1, D-86135 Augsburg, Germany 
(Dated: February 6, 2008) 

Abstract 

We investigate a subdiffusive, fractional Fokker-Planck dynamics occurring in time-varying po- 
tential landscapes and thereby disclose the failure of the fractional Fokker-Planck equation (FFPE) 
in its commonly used form when generalized in an ad hoc manner to time-dependent forces. A 
modified FFPE (MFFPE) is rigorously derived, being valid for a family of dichotomously alter- 
nating force-fields. This MFFPE is numerically validated for a rectangular time-dependent force 
with zero average bias. For this case subdiffusion is shown to become enhanced as compared to 
the force free case. We question, however, the existence of any physically valid FFPE for arbitrary 
varying time-dependent fields that differ from this dichotomous varying family. 

PACS numbers: 05.40.-a, 05.40.Fb, 05.60.-k, 02.50.Ey 
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Normal Brownian motion occurring on potential landscapes that vary in time is known to 
exhibit a multifaceted collection of interesting phenomena, such as Brownian motors, anoma- 
lous nonlinear response behaviors, and stochastic resonance [l], to name a few. Therefore, 
it is tempting to ask, whether an explicit time-dependent force entails a similarly versatile 
scenario also in the case of anomalously slow relaxation processes, relevant in many systems, 
such as polymer chains, networks, proteins, glasses and charge- carriers in semiconductors 
Q]. This issue is in fact contained already in the first works on the motion of charge- carriers 
in semiconductors 3|, and has been the subject of some further investigations ever since, see 
e.g. the works {4], 5 Ja, Q], but never really has attracted proper attention on its fundamental 
level. Ultraslow relaxation in time- dependent external potential-fields thus still constitutes 
a challenge that is far from trivial. 

A widely used approachto study subdiffusive processes is based on the fractional Fokker- 
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Here F(x) is the force, rj a is the fractional friction coefficient, K a is the fractional free diffusion 
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coefficient, and oD t denotes the Riemann-Liouville fractional derivative, 
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For time-independent forces the FFPE ([I]) can be rigorously derived from continuous time 



random walk (CTRW) theory 



, |9] , which corresponds to the simple random walk, including 



the non-Markovian character of the process via a Mittag-Leffier residence time distribution 
(RTD) V>(t) oc r~ x ~ a (0 < a < 1) [a,la,llQ|. 



With this work we show that the FFPE (CQ) fails in the case of a time- dependent force 
F(x,t). Furthermore, we argue that the FFPE ([I]), when generalized ad hoc to a time- 
dependent force, does not correspond to a physical stochastic process. This affects a still 
steadily growing body of current research [11], and implies that these so obtained results 
therein are physically defeasible. 

In different context, the study of a subdiffusive dynamics in the case of a purely time- 
dependent force F(t) has given rise to a fractional Fokker-Planck equation which differs from 
Eq. ([1]) 7]. In this letter, we derive an equation of similar form for the class of dichotomously 
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alternating force-fields F(x,t) = F(x)£(t) with = ±1, varying in space and time. In the 
case of a Mittag-Lefner RTD it reads, 
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Below, we prove this form in terms of CTRW theory and additionally validate its correctness 
via the comparison of the analytical solutions of this so modified fractional Fokker-Planck 
equation (MFFPE) (j3J) for a rectangular time- varying periodic force F(x, t) = F(t) = £(t)F 
with the numerical simulations of the underlying CTRW. Different force-fields, such as 
F(x, t) = ± sin(x) have also been successfully tested (the details will be presented in a 
longer follow up work). Our main point is, however, that the reasoning provided in proving 
([3]) forces us to scrutinize the physical validity of this so modified FFPE in ([3]) already be- 
yond a dichotomous driving F(t); as e.g. it happens for a sinusoidal driving £(£) as used in 
Ref. fl. 

As an interesting result, we also show that a symmetric dichotomous force with average 
zero bias enhances the diffusion in respect to the free case. Furthermore, it is found that for 
sufficiently slow driving the effective fractional diffusion coefficient exhibits a maximum 
vs. the fractional exponent a. 

It is well known that neither a non-Markovian Fokker-Planck equation nor its solution 
with the initial condition P(x,t) = 5(x — Xq) (i.e. the two-event conditional probability 



P(x,t\xo)) can fully define the non-Markovian stochastic process 12j. This is due to the 
fact that all non-Markovian processes, such as a CTRW with a non-exponential waiting time 
distribution, lack per definition the factorization property, which would allow to express all 
the higher-order (multi-event) probability density functions in terms of the first two ones. 
Because a CTRW is at the root of the FFPE (CD) 9], in order to generalize the latter to 
the time-dependent forces, one again starts from CTRW theory. However, the usual scheme 
of merely replacing a time-independent force F in Eq. ([1]) in an ad hoc manner with a time- 
dependent F(t) is doomed to failure. The reason is that the underlying CTRW possesses 
a RTD with an infinite mean. Thus, any regular driving with a large but finite period is 
nonadiabatic. This very circumstance lies at the heart of the overall failure of Eq. (OQ) for 
time-dependent force fields. 

In terms of a renewal description, a CTRW is a semi-Markovian process, meaning that 
the sojourn times spent on the localization sites are independently distributed. Let us 



consider a one- dimensional CTRW on a lattice Xi = iAx (i = 0, ±1, ±2, . . . ). After a time 
r drawn from the RTD ipi(r), the particle at site i jumps with the probability qf to one 
of the nearest neighbor sites. The external force-field F(x) specifies both i/j^t) and qf, see 
Ref. 9]. Modulating the force F(x) in time, qf assume obviously a time-dependence and 
ijjf(T\t) = qf(t + r)^j(t + r, t) become conditioned on the entrance time t for the site i jl3 |. 
For a Markovian CTRW with time- dependent rates wf(t) it is known that 
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ipi(t + r,t) = Wi(t + r) exp 



Wi(t')dt' 
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with Wi(t) = wf(t) + w^(t) and qf(t) = wf(t)/wi(t). For a driven non-Markovian CTRW, 
however, a relation similar to Eq. (j3J) is lacking. As a result, the use of a FFPE when 
generalized to the time-dependent case of a time- varying force-field remains moot. The 
usual scheme of the derivation of the generalized FFPE from the underlying CTRW can be 
used only if ipi{r) remains unmodified by the time-dependent fields, i.e. if only the jump 
probabilities qf (t) change. We consequently find that ipf(r\t) = qf(t + r)^j(r). Thus, the 
RTD ipi(r) remains unaffected only in the case of a dichotomous flashing force F(x,t) = 
F(x)l;(t), where = ±1 is a general dichotomic function of time t which can change 
periodically or also stochastically. Then, qf(t) = exp(F(xi)t;(t)Ax/2)/[exp(F(xi)Ax/2) + 
exp(— F(xi) Ax/2)]. We assume that F(x) is continuous. Then, the MFFPE ([3]) can be 
derived rigorously in the continuous space limit. The derivation precisely follows the same 
reasoning as detailed in Ref. [3], while taking ipi(r) as being the Mittag-Lefner distribution. It 
must be emphasized that for other driving forms £(£), e.g. for a sinusoidal driving Fq sin(u)t), 
this outlined derivation becomes flawed because ipi(r) is affected by such time- varying fields, 
as being unveiled already with Eq. (Tj0). We remark also that due to the weak ergodicity 
breaking [jj], [jjj] this MFFPE (J3j) describes the dynamics of an ensemble of particles rather 
than the dynamics of an individual particle. 

Notably, Eq. (JTJ applied to the case of a time-dependent force may well define an interest- 
ing mathematical object in its own right; its connection to a known physical process, however, 
remains open to question. One may attempt to justify Eq. ([T]) for a time-inhomogeneous 
situation by appealing to the concept of " subordination" ; noting that this equation corre- 
sponds to a random process described by a usual Langevin equation but with the operational 



time being a random stochastic process 16|, A time-dependent physical force, however, 



varies in deterministic, real time, which physically cannot be transformed to random time. 



In the following we study the particular case of a periodic rectangular driving force 
F(t) = Fo(— iy 2t / T0 \ where To denotes the time-period and [a] is the integer part of a. Put 
differently, we consider a dichotomous modulation of a biased free subdiffusion where the 
absolute value of the bias is fixed, but the direction of the force flips periodically in time. 
The average bias is zero. Let us begin by finding the recurrence relation for the moments 
(x n (t)). Multiplying both sides of Eq. (J3j) by x n and integrating over the x-coordinate one 
obtains, 

= nFm-^it))/^ 

+ n{n-l)K a D\- a {x n - 2 {t)), (5) 

for > 2. When n — 1 the last term on the right hand side of the latter equation is absent. 
Then, 



d(x(t)) = F(t) Dl _ ai = F(t) e _ x _ 
dt rj a 1 VaF(a) 

Integrating Eq. in time, the solution for (x(t)) reads: 
(x(t)) 
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Here, v a = F /rj a and counts the number of time periods passed. The analytical solu- 
tion ([7]) for the mean particle position (x(t)) from the MFFPE ([3]) is compared with the 
numerical solution of the CTRW in Fig. [1] for different values of the fractional exponent 
a. The good agreement between our analytical and numerical results confirms that Eq. (J3]) 
is a correct method to describe the CTRW driven by a rectangular time- dependent force. 
Furthermore, the results depicted in Fig. [T] exhibit the phenomenon of the "death of linear 
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FIG. 1: (Color online) Average particle position (x(t)) for various values of the fractional expo- 
nent a: Symbols represent the numerical results for the CTRW obtained by averaging over 10 6 
trajectories, while continuous lines represent the analytical solution (|7|) of the MFFPE ([3]). The 
time-period of the force is to = 1 and Fn[(r] a J~K^i) = 1 is used in numerical simulations. The 
simulation algorithm is described in Ref. [if 



response" of the fractional kinetics to time-dependent fields in the limit t — > oo, reported also 
in Refs. {5, 7|; i.e. in the long-time limit the mean particle position approaches a constant 
value, rather than being oscillatory, i.e., 



(ar(oo)) = v a T«b{a)/T{a + 1) 



(10) 



where b(a) = ^^LoP( n + 1/2)" — n a — (n + l) a ], with the amplitude of the oscillations 
decaying to zero as l/t 1-a , see Eq. (jSJ). The function b(a) describes the initial field phase 
effect which the system remembers forever when a < 1. It changes monotonously from 
6(0) = 1 to b(l) = 0. The averaged traveled distance (x(oo)) scales as Tq = (2n/Q) a , where 
Q is the corresponding angular frequency. This "death of linear response" to time-periodic 
fields is also in agreement with the results for a driven non-Markovian two state system [l^] 
in the formal limit of infinite mean residence times. 

We next study the mean square displacement and the effective fractional diffusion co- 
efficient ■ We recall that the free fractional diffusion is described by (5x 2 (t)) = 
2K nt a / T(l + a) oc t a , while in the presence of a constant bias, surprisingly, (5x 2 (t)) oc t 2a 
[3I. Tlo|| . Strikingly enough, the same difference in the behavior remains true for a CTRW 



proceeding in a periodic potential with zero bias for which (Sx 2 (t)) oc t a 19( and in a wash- 
board potential with finite bias for which (5x 2 (t)} oc t 2a [^[ijj], described by the FFPE. The 
question thus arises, whether a time-modulated subdiffusion follows the biased fractional 




FIG. 2: (Color online) Scaled effective fractional diffusion coefficient Ka versus fractional ex- 
ponent a for different driving periods tq. The analytical prediction (I14p (continuous lines) is 
compared with the numerical results (symbols) obtained from the CTRW by averaging over 10 5 
trajectories. For tq > 2ir exp[— 4- In <?(a)| a =o] ~ 8.818 the effective fractional diffusion coefficient 



(eff) 



(a) exhibits a maximum. 



subdiffusion behavior oc t 2a , or rather assumes the unbiased behavior oc t a , as the average 
bias is zero. To obtain the answer we use the Laplace-transform method and the Fourier 
series expansion for the driving force F(t + r ) = F(t) with frequency Q = 2tt/t , 



F(t) = J2 fn****, f- n = r n 



(11) 



with / 2n = and /2n+i = — (2z/7r)Fo/(2n+l) for the rectangular driving under consideration. 
We assume that (x(0)) = and (x 2 (0)) = and denote the Laplace-transforms of the first 
and second moment by x(s) and y(s), yielding 



sx 



(12) 



where v n = f n /v&- Because (x(oo)) is finite, we evaluate the asymptotical behavior of (x 2 (t)) 
rather than (8x 2 (t)). Because the Laplace-transform of Dl~ a (x(t)) is s~ a sx(s), the Laplace 
transform of the second moment reads 



sy(s) = 2n a /s a + 2 v m (s — imfl)' 

m=—oo 

oo 
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Note that in the double sum only the terms with m = —n contribute to the effective 
subdiffusion coefficient, being averaged over the driving period. Therefore, we find that 



is a function decaying from g(0) = 1 towards g(l) = and ((x) is the Riemann's zeta- 
function. From Eq. f|T3|) one finds that the asymptotic behavior of the mean square dis- 
placement is proportional to t a as in the force free case. It is now characterized, how- 
ever, by an effective fractional diffusion coefficient Ka S ^ instead of the free value K a , i.e., 
(5x 2 (t)} = 2Ka S ^t a /r(l + a) for t — > oo. The driving-induced part of the effective subdif- 
fusion coefficient is directly proportional to the square of driving amplitude and inversely 
proportional to Q a . For slowly oscillating force fields this leads to a profound acceleration 
of subdiffusion as compared with the force free case, see in Fig. [2j An optimal value of the 
fractional exponent a exists, at which the driving-induced part of the effective fractional 
diffusion coefficient possesses a maximum. 

In this letter we discussed the dynamics of anomalously slow processes in time-varying 
potential landscapes within the CTRW and FFPE descriptions. We demonstrated that the 
common form of the FFPE given by Eq. ([I]) is not valid for time-dependent forces; it fails to 
correspond to the underlying CTRW modulated by an external time-dependent force field. 
A modified form of the FFPE, Eq. ([3]), is derived for dichotomously alternating force-fields. 
As an exactly solvable example we studied a periodic rectangular force with zero average 
and successfully tested the analytical results via numerical simulations of the underlying 
time-modulated CTRW. 

Our study, however, is not able to validate the correctness of the MFFPE fl3]) when 
extended ad hoc to an arbitrary time-dependent potential landscape different from the di- 
chotomous case. A description of time-dependent fields via subordination in conjunction 
with a CTRW approach is also doomed to failure because of the distinct difference between 
the deterministic physical time and the merely mathematical random subordination time. 
As a matter of fact, any slowly non-zero frequency time- varying force varies infinitely fast 
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where 



g(a) = (2/tt 2 )C(2 + a) [4 - 2~ a ] cos(7ra/2) 



(15) 
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within the realm of fictitious, operational subordination time. This causes CTRW subdif- 
fusion to fail in responding to time-periodic fields in an ordinary manner. In addition, all 
those theories modeling dielectric response which are based on such an approach are thus 
also physically defeasible. A way out of this dilemma consists in relying on models of driven 
sub diffusion which either are based on the generalized Langevin dynamics [20(] or on frac- 
tal Brownian motion. The challenge of modeling subdiffusion in a time-varying potential 
landscape thus necessitates plenty of further enlightening research. 
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